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SUMMARY 

The effects of existing frictional heating were analyzed to deter- 
mine the conditions under which ice formations on aircraft surfaces can 
be prevented. A method is presented for rapidly determining by means 
of charts the combination of-Mach number, altitude, and stream tempera- 
ture which will maintain an ice-free surface in an icing cloud. The 
method can be applied to both subsonic and supersonic flow. The charts 
presented are for Mach n\imbers up to 1.8 and pressure altitudes from sea 
level to 45,000 feet. 


INTRODUCTION 

As the operational speed of aircraft is increased through the tran- 
sonic region, the frictional heating available to prevent the formation 
of ice on the aircraft becomes an Important quantity. The operation of 
missiles and interceptor aircraft is therefore possible under pre- 
determined values of altitude, stream temperature, and Mach number which 
will preclude icing. The set of flight conditions which provides a 
sxirface temperat\ire of 32° F for a particular point on a body traveling 
in an icing cloud is termed the "icing limit" for that point. The 
analysis presented in reference 1 relates the frictional heating, evapor- 
ation, and heat transfer for this particular condition. This analysis 
was applied to a symmetrical diamond profile airfoil at zero angle of 
attack and it was found that the most critical region of the airfoil in 
terms of icing limit was on the rear surface behind the midchord 
(shoulder). Although this region is not subject to direct impingement 
of water droplets, the water film caused by droplets impinging on the 
front s\irface runs, back over the shoulder and extends the hazard of 
icing over a large percentage of the airfoil. Icing might therefore 
occur at any point on a body which is subject to droplet impingement or 
njnback. 

Experimental verification of the calculations presented in refer- 
ence 1 was obtained at the NACA Lewis laboratory on a diamond airfoil of 
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5 50 inch chord and 6.82 percent thickness for a Mach numiber of 1.36 and 
several pressure altitudes . The good agreement between experimental and 
analytical results provides credence for the analytical method presented 

therein. 

The procedure used for calculating the icing limit in reference 1 
can be applied generally to obtain the icing limit for any body provided 
the flow conditions about the body are known. The charts presented 
herein are based on the analysis of reference 1 for calculating the icing 
limit at any point on a body. These charts are similar to that shown in 
figure 1 of reference 1 but cover an altitude range from sea level to 
45 000 feet . Several auxiliary charts are also presented from which the 
pressure coefficient, velocity ratio, or local Mach nuaiber at a particu- 
lar point on the body can be converted to the parameters utilized in the 
icing limit charts. 


^ ■ l^^OLS 


The following symbols are \ised in this report : 

Pi - Po 


pressu-ve coef ficien'b y 


1/2 PoV 



L 

M 


“a 

“e 

P 


r 

T 


specific heat of air at constant pressure, Btu/(lb)(°F) 

vapor pressure, Ib/sq ft 

coefficient of evaporation 

coefficient of heat transfer 

latent heat of vaporization, Btu/lb 

Mach nimiber 

molecular weight of air 
molecular weight of water vapor 
static pressure, Ib/sq ft 
recovery factor 
static temperature, °R 


OU minimum free- stream static temperature corresponding to ice-free 
condition on surface as defined in eq. (2) , R 
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V velocity, ft /sec 

r ratio of specific heats of air, 1.400 

p density of air, slugs /cu ft 

Subscripts : 

0. free-stream static conditions 

1 local conditions at edge of boundary layer 

s surface 

ANALYSIS 

The generalized heat balance in an icing cloud for an insulated 
body heated only by frictional effects is given as follows : 

(l) Heat due to the frictional, or viscous, effects 

plus 

(2) Heat due to the kinetic energy of the water droplets 

plus 

(3) Heat of fusion 
equal 

(4) Heat lost by convection 
plus 

(5) Heat for evaporation of water 
plus 

(6) Heat required to raise temperature of water droplets from stream 

temperature to surface temperature. 

Hardy's relation (ref. 2) is obtained by equating the frictional 
term (l) to the sum of the convective term (4) and the evaporative term 
(5). The heat of fusion term (3) is equal to zero if an ice-free 
sijrface is defined as being fully wet at 32° F with no Ice particles in 
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"the surface water film* In addition, terms (2) and (6) are of the same 
order of magnitude and tend to offset each other for calculations of the 
icing limit. Both increase in value with increasing tfach number; large 
values of either term are found only at large values of liquid-water 
content and flight speed. Such a combination is highly unlikely in 
actual flight since high speeds usually occur at high altitudes, and 
high altitudes imply small liquid-water contents with the exception of 
flight in cumulus clouds of great vertical extent. 

The results of reference 1 indicate that the relation of Hardy 
(ref. 2) applies in particular to a surface barely wetted by a very thin 
film, but is probably nearly correct for the whole range of speeds and 
liquid-water contents of practical interest for high-speed flight. 

The relation of Hardy is given in reference 1 as 


n . r(r-l) 2 

ke“eL / eg 

1 + 2 % 

' ^h®a°p W ■ °s Pi - °1/ 


Because interest in this problem is fixed on a definition^of flight cir- 
cumstances that provides local s\irface temperatures of 32 F, the terms 
T , L, and e are constants. Hardy observed that for the range of 

temperatures near 32° F the ratio of the evaporation coefficient kg to 
the heat-transfer coefficient kj^ is very nearly unity. ' In addition, 

the air in a cloud is fully saturated at the static or free-stream condi- 
tions and the free-stream vapor pressure eg is therefore equal to the 

saturated vapor pressure at the free-stream static temperature . If it 
is assumed that the flow about the body, outside the boundary layer, is 
accomplished with no change in phase, that is, no condensation or evap- 
oration, then Dalton's law of partial pressures applies and 

®0 Po 


Equation (l) can therefore be rewritten in the form 


Tl( 

1 + r 

12.75 

°0 


in n 

m 0,0 

^0,c 

(Pi/Pq) Pq ” 12.75 

Pq ■ °0_ 


( 2 ) 


which is quite convenient for purposes of calculation because 
^i/Tq ^ 1 + r ] and Pi/pq can be readily calculated at any 
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point on a body for shock-free flow from the pressure or velocity dis- 
tribution, the stream Mach number, and the recovery factor r. When the 
flow about the body contains a shock, then the parameters 


('TiAo,c) ^ Pl/Po can be obtained from the local 

Mach number, the stream Mach number, the pressure coefficient, and the 
recovery factor r. 


PREPARATION OF ICING LIMIT AND AUXILIARY CHARTS 

Because Bq is a function of the solution of equation (2) 

must be accomplished by trial and error. In an effort to reduce the 
laboriousness of such calculations a number of charts have been made and 
are presented in figure 1. For each pressure altitude considered, the 
free-stream static temperature corresponding to an ice-free condition 

Tq,c plotted as a function of the parameter (^±/^q ^1 + r Mi2^ 

by means of equation (2) for constant values of P^^/Pq* For each set of 

values of Pj^/Pq ^ "*■ ^ ^ ^^®®“Stream static 

tenqperatvire exists which corresponds to a surface temperature of 32° F 
and therefore represents the minimum free-stream static temperature for 
an ice-free surface The lower limit of the free-stream static 

temperature was considered to be -40° F, since it has been shown (refs. 3 
and 4) that supercooled water droplets are not likely to exist below this 
temperatxire . 


The ten icing limit charts presented (figs. l(a) to l(j)) were cal- 
culated for pressure altitudes from sea level (fig. 1(a)) to 45,000 feet 
(fig. l(j)) in 5000 foot increments. The use of the icing limit charts 

requires the knowledge of the parameters ^ r 

Pl/po* For shock-free flow both these parameters can be readily calcu- 
lated from the pressure coefficient Cp or the ratio of local to free- 
stream velocity Vq^/Vq and the stream Mach number. For most subsonic 


airfoils either the experimental pressure coefficient or the calculated 
velocity ratio is known. A variety of methods is available in the lit- 
erature for determining local values of Cp or Vj^/Vq for almost any 

body in shock-free flow (refs. 5 to 9) . 


In the transonic-supersonic regime, wherein a shock normally exists 
in the flow field, the flow field about a body is not defined by the 
pressure coefficient and stream Mach number alone and an additional 
parameter such as shock wave angle or local ^fe.ch number is required. 



6 


NACA TN 2914 


Most of the results in the literature are presented in terms of local 
pressure coefficient and local Mach number for constant values of free- 
stream Mach number (ref . lO) • 


The pressvire coefficient may be related to the pressure ratio and 
the stream Mach number by the compressible flow relations for a perfect 
gas without the assumption of isentropic flow as follows: 


1 


Pi - Po 2(pi - Pq) 


•* 0 V 
2 ^0 0 


r Po Mq 






^ = 1 + I Mq 2 Cp = 1 + 0.7 Mo’-^^ Cp 


(3) 


Hence this result applies to both shocked and shock-free flows. 


The velocity ratio may be related to the pressixre ratio and the 
stream Mach number in the following manner for shock-free flow: 



PlVi Po 




pq Pi 
^1 ^0 ^0 



r-i 

r 



The static pressure - Mach ntunber relation for isentropic flow is given 



1 

1 + 


Solving for 



gives 
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*1 2.236 

Vo " Mo 


1 + 0.2 Mq 


Pq. 


0.286 


(4) 


The parameter r is related to P^/Pq^ “the 

stream Mach number, and the recovery factor in the following way for 
shock-free flow; 


T. 


T. 


0,c 


1 + ^ r M- 


r-1 

r 


.1-© (■•¥ 


r Mn 


hut 


r-1 


y-1 2 

1 + V“i 


■© (‘*¥^1 


Multiplying hoth sides hy r, adding unity to each side, and rearranging 
terms yield 


r-1 

r 


1 + r 


(l + Iji Mo=) r + (1 - r) 


Therefore, hy substitution. 


( 

= (1 + 0.2 Mq^) r + (1 - r) 


Hi 

r 


0.286 


(5) 



8 


MCA TN 2914 


For shocked flows the pressure ratio can be obtained by equation (3)^ 

but the parameter (^1 + r can no longer be related to 

the pressure ratio and stream Mach nxanber as for shock-free flows 
(eq. (5)). 

The parameter ^Ti/Tq^c) ^ *^1^) however, be related 

quite simply to the local and free-stream Mach numbers. The temperatiore 
ratio 1-^/Iq is given in terms of Mach number by the energy relations 

which apply to both shocked and shock-free flows as follows: 


1 

1 + ^ Mi^ 


therefore 


T 


1 




1 + ^ Mo^ 

1 + Mi^ 



( 6 ) 


The pressure ratio P]_/pq can be readily obtained by use of equa- 
tion (3) if Mq and Cp are known. The term Mq Cp can be calculated 

from experimental pressure distributions which are usually presented for 
each value of Mq in terms of Cp. 


Theoretical velocity distributions for many bodies are given in terms 
of the velocity ratio or the square of the velocity ratio. The conversion 
from velocity ratio to pressure ratio (eq. (4)) is shown in figure 2, 
where the pressure ratio is plotted as a function of velocity ratio for 
constant values of the stream Mach number. 


Figure 3 shows the conversion from pressure ratio to the parameter 
^i/Tq ^ ^1 + r Mq^^ as given by equation (5) for shock-free flows. 

The parameter ^^l/^O c^ ^2 ~ ^ plotted as a function of the 

stream Mach number Mq for constant values of the pressure ratio V-^/Vq’ 

Figures 3(a), 3(b), and 3(c) are given for values of the recovery factor 
r of 0.80, 0.85, and 0.90, respectively. 


The relation given by equation 
the parameter (^i/Tq,c^ ^ 


(6) is shown plotted in figure 4 where 
is plotted as a function of the 
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stream Mach mmiber for constant values of the local to free- stream Mach 
number ratio M^^/Mq. Figures 4(a) ^ 4(h), and 4(c) are given for values 

of the recovery factor of 0-80, 0-85, and 0-90, respectively. 


Discussion and Illustrative Examples 

Most aircraft components which are exposed to icing are of airfoil 
or streamlined shape. For this general class of bodies most of the 
cloud water droplets impinge near the leading edge or nose region. The 
icing hazard, however, is not limited to the region of direct impinge- 
ment, since if the impinged cloud droplets do not freeze (because of 
frictional heating) they create a surface water film which flows back 
from the impingement area. The whole region wetted by the film must 
therefore be considered as a region of possible icing hazard. 

The critical point on an airfoil or similar body (i.e., the point 
where the largest value of Tq ^ is required for a particular value of 

stream Mach nimiber) occurs at the minimum pressure point because the 
decreased pressure enhances the evaporation and results in a reduced 
surface temperature. If impingement occurs only near the leading edge 
or nose region, the evaporation in the region ahead of the point of 
minimimi pressure may be sufficient to remove all the water film and, 
in such a case, the results obtained from charts l(a) to l(j) would not 
s-PPly- The actual critical point might therefore occur foi*ward of the 
minimum pressiire point; therefore^ values of Tq^^ based on minim\im pres 

sure at a particular Mach mmiber would be greater than actually required. 

, The use of charts 1(a) to l(j) can be shown most readily by means of 
several examples. The following example illustrates typical results in 
the subsonic speed range. 

Example I . - Calculation is made of the free-stream static temper- 
ature required for an ice-free surface as a function of Mach number for 
an NACA 65-206 airfoil at 15,000 feet pressure altitude and 1^ angle of 
attack - 

The velocity ratio V^^/Vq can be found for this airfoil by the 

methods and results presented in reference 5. The maximimi velocity ratio 

(minimum pressure ratio) for this airfoil was determined to be 

V^/Vq = 1.139 and occurs at the 45-percent chord station on the upper ' 

surface. Assuming a value of the recovery factor of 0.85 and that the 
surface is wet at this point determines Tq^^ in the following way: 
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Free -stream 
Mach number, 

«0 

Static-pressure 
ratio, P-j^/Pq 

(fig. 2(b) and 
Vi/Vq) 

Parameter 

To,c(' ) 

[fig. 3(b), M^ 

Minimum free-stream 
static temperature, 

'^0,c^ 

OR 

[fig. 1(d)] 

0.400 

0.970 

1.026 

486.8 

.500 

.950 

1.041 

483.8 

.600 

.929 

1.059 

480.0 

.700 

.902 

1.079 

475.6 


The values of Tq ^ obviously decrease with increasing Mach nuniber as 

wo\ild be expected. The results of example 1 show, however, that even 
at a Mach number of 0.7, protection is provided at stream temperatures 
greater than 16° F. 

Example 2. - This example is presented to show typical calculations 
and results in the transonic speed range. Determination is made of the 
relation between stream Mach nturiber and stream static temperature which 
will provide an ice-free surface at the midchord of an 8.8 percent thick 
circtaar arc airfoil at zero angle of attack for altitudes of 10,000, 
25,000, and 40,000 feet. The recovery factor r is 0.90. 


The values of pressiare coefficient Cp and local Mach nxmiber Mp 

can be obtained from reference 10 for stream Mach numbers of 0.848 to 
1.500 and are listed as follows: 


Mo 

Ml 


0.7 

Pi 


^ (A)l 

/0,c 

[fig. 4(c), 
\ MpX 

T0,c. 

(ref. 

10) 

(ref. 

10) 

at 

mid- 

chord 

(ref. 

10) 

at 

mid- 

chord 

Cp 

(cal- 

cu- 

lated) 

PO 

Mo 

Altitude, 

ft 

(eq. 

(3)) 

(cal- 

cu- 

lated) 

10,000 

(fig- 

1(c)) 

25,000 

(fig. 

1(f)) 

40,000 

(fig. 

l(i)) 

0.848 

1.023 

-0.355 

-0.179 

0.820 

1.206 

1.123 

460.0 

472.5 

484.0 

.935 

1.140 

-.330 

-.202 

.798 

1.219 

1.150 

4:52.0 

466.5 

481.5 

1.110 

1.160 

-.045 

-.039 

.961 

1.045 

1.220 

423.5 

440.0 

464.0 

1.160 

1.183 

-.028 

-.026 

.976 

1.020 

1.241 


432.5 

459.0 

1.200 

1.195 

.010 

.010 

.990 

.996 

1.260 


425.5 

454.0 

1.250 

1.350 

1.225 

1.306 

.025 

.028 

.027 

.036 

1.027 

1.036 

.980 

.967 

1.282 

1.331 



447.5 

435.0 


^A = 1 + r . 


i 
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The results of example 2 are shown plotted in figure 5, where the 
stream static temperature corresponding to an ice-free condition is 
plotted as a function of the fli^t Mach number for the three altitudes 
considered. The effect of altitude is clearly illustrated. The free- 
stream static temperatures required for protection axe higher at the 
higher pressure altitudes, as would be expected; since, if other condi- 
tions are equal, the decreased pressure enhances the evaporation rate. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 4, 1952 


BEFERENCES 

1. Callahan, Edmimd E., and Serafini, John S. : Analytical Investigation 

of Icing Limit for Diamond Airfoil in Transonic and Supersonic Flow. 
NACA TN 2861, 1953. 

2. Hardy, J. K. : Kinetic Temperature of Wet Surfaces - A Method of Cal- 

culating the Mount of Alcohol Required to Prevent Ice, and the 
Derivation of the Psychrometric Equation. NACA ARR 5G13, 1945. 

3. Schaefer, Vincent J.: The Production of Clouds Containing Supercooled 

Water Droplets or Ice Crystals Under Laboratory Conditions. Bull. 
M. Meteorological Soc., vol. 29, no. 4, April 1948, pp. 175-182. 

4. Dorsch, Robert G., and Hacker, Paul T. : Photomicrographic Investiga- 

tion of Spontaneous Freezing Temperatures of Supercooled Water 
Droplets. NACA TN 2142, 1950. 

5. Abbott, Ira H., von Doenhoff, Albert E., and Stivers, Louis S. , Jr.: 

Summary of Airfoil Data. NACA ACR L5C05, 1945. 

6. Theodorsen, T., and Garrick, I. E. : General Potential Theory of 

Arbitrary Wing Sections. NACA Rep. 452, 1933. 

7. Allen, H. Julian: A Simplified Method for the Calc^llation of Airfoil 

Pressure Distribution. NACA TN 708, 1939. 

8. Jones, Robert T., and Cohen, Doris: A Graphical Method of Determining 

Pressure Distribution in Two-Dimensional Flow. NACA Rep. 722, 1941. 

9. Kaplan, Carl: Potential Flow about Elongated Bodies of Revolution. 

NACA Rep. 516, 1935. 

10. Bryson, Arthur Earl, Jr. : An Experimental Investigation of Transonic 

Flow Past Two-Dimensional Wedge and Circular-Arc Sections Using a 
Mach-Zehnder Interferometer. NACA TN 2560, 1951. 



12 


NACA TN 2914 



Figure 1. - Variation of minimum free- stream static temperature Tq ^ corresponding 
to the Ice-free condition with parameter (T^/Tq q)(^ ^ various 

pressure ratios (eq. (2)). 
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Figure 1. - Continued. Variation of minimum free-stream static temperature T* cor- 

u, c 

responding to the ice-free condition with parameter (Tq^/Tq^^) r for 

various pressure ratios (eq. (2)). 



Mlninium free-stream static temperature. 
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480 


*1 470 


460 



450 


1.32 


Parameter (T^/Tq^c^ ''' ^ ^1^) 

(c) Pressure altitude, 10,000 feet. 

Figure 1. - Continued. Variation of minimum free-stream static temperature corresponding to 

the ice-free condition with parameter (T^/Tq q) ^ various pressure ratios 

(eq. (2)). 





free- stream static temperat\xre 
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(d) Pressure altitude, 15,CXXD feet. 

Figure 1. - Continued. Variation of minimum free-stream static temperature corresponding to 

u,c 

the Ice-free condition with parameter (T^/Tq q) ^ various pressure ratios 

(eq. (2)). 
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Figure 1. - Continued. Variation of m i ni mum free-streeun static ten5)erature Tq ^ corresponding to the ice-free condition 
with parameter (T^/Tq^^)(i + ^ r ^or various pressure ratios (eq. (2))! 
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Figure 1. - Continued. Variation of minlnium free-stream static temperature Tq^^, corresponding to the ice-free condition 
with parameter (T^/Tq q) ^ various pressure ratios (eq. (2)). 
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(h) Concluded. Pressure altitude, 35,000 feet. 

Figure 1. - Continued. Variation of minimum free-stream 
static temperature Tq - corresponding to the ice -free 

condition with parameter (^i/Tq^q) r j for 

various pressure ratios (eq. (2)). 
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Figure 1. - Continued. Variation of minimum free-stream static temperature Tq ^ corresponding to the Ice-free condition 
vlth peirameter (Ti/Tq „) r for various pressure ratios (eq.. (2)). 
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(i) Concluded. Pressure altitude, 40,000 feet. 


Figure 1. - Continued. Variation of Tnlnlmnm free-stream static 
temperature corresponding to the Ice-free condition with 

parameter (Ti/To,c) ^ for various pressure ratios 

(eqi. (2) ) . 
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Minimum free-stream static teBg[)erature^ 
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(j) Concluded. Pressure altitude, 45,000 feet. 

Figure 1. - Concluded. Variation of minimum free-stream static 
temperature corresponding to the ice-free condition 

with parameter ^1 + r for various pressure 

ratios (eq.. (2)). 
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Figure 2. - Variation of static -pres sure ratio with velocity ratio in shock- free flow for various free-stream Mach numbers 
(eq. (4)). 
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fb) Velocity ratios greater than 1. 

Figure 2. - Concluded. Variation of static -pressure ratio with 
velocity ratio in shock-free flow for various free-stream Mach 
numbers (eq. (4)). 
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Figure 5. - Variation of minimum free-stream static temperature for ice-free surface 
as function of free-stream Mach number for station at 60 percent chord of 8.8 percent 
thick circular arc airfoil. 
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